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We study weakly coupled antiferromagnetic spin chains in their ordered phase by combinining an 
exact solution of the single-chain problem with an RPA analysis of the interchain interaction. A 
single chain is described by a quantum Sine-Gordon model and dynamical staggered susceptibilities 
are determined by employing the formfactor approach to quantum correlation functions. We consider 
both antiferromagnetic order encountered in quasi-lD materials like KCUF3 and spin-Peierls order 
as found in CuGeOa. 
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I. INTRODUCTION 

The increasing amount of neutron data on quasi-lD antiferromagnets calls for creation of a theory capable of 
describing the behaviour of such strongly anisotropic nonlinear systems in greater detail. Of particular interest with 
regard to recent experiments on CuGeOa and KCUF3 is the structure of the multiparticle continua. In the present 
paper we discuss correlation functions of three closely related problems involving spin-i Heisenberg chains: (i) a single 
chain with a static alternating exchange, (ii) a single chain in a staggered magnetic field, (iii) a system of weakly 
coupled spin-i Heisenberg chains at low temperatures in the ordered phase. In the latter case we consider two types 
of order: (a) spin-Peierls 0-^] and (b) antiferromagnetic . 

An experimental realisation of situation (i) may be achieved in (VO)2P20y. The analysis of the recent experiments 
conducted by indicates that the strongest exchange occurs through alternatively arranged molecules of two different 
types such that the ratio of the exchange integrals is J'/ J « 0.722. 

The second situation (ii) may be realized when a Heisenberg magnet with a staggered Lande factor is placed in an 
external magnetic field. This is the case for copper benzoate |9[| but the difference in Lande factors is very small and 
for reasons that will be discussed elsewhere the theory presented here does not apply. 

The spin-Peierls transition is a magnetoelastic transition which occurs in quasi-lD antiferromagnets due to an 
effective four-spin interaction. Such interactions may be generated by phonons which modify the exchange integrals; 
however, the virtual absence of softening of the phonon spectrum in the spin-Peierls material CuGe03 JlOj ] suggests 
that there may be also other mechanisms. The interaction between the staggered parts of energy densities is strongly 
relevant and resolves in dimerization of the lattice at a certain temperature T c and the formation of spectral gaps 
in the spectrum of magnetic excitations. In any realistic system phonons are three-dimensional which determines a 
three dimensional nature of the spin-Peierls transition. A simple model with a phonon mechanism has the following 
Hamiltonian 



H = (1 + (-l)"u„,r) S„, r ' S n+ i, r + i u n,vD{v - r')u n , T > 

r,n n,rr' 

+ a J ^ n ' r ' S"+ 2 . r + 51 K(c) S n ,r • S„, r +c , (1) 



r,n r,c,n 



where n and r label lattice sites along and perpendicular to the chains and c are vectors connecting neighbouring 
chains. In order to make contact with experiments we have included n.n.n. intrachain as well as interchain magnetic 
interactions. 

Since J is usually much smaller than the phonon frequency, one can neglect the kinetic energy of the w-field and 
consider u as commuting numbers. One can denote 

e(n,r) = (-l)"(S„, r ■ S„+i, r ) (2) 
and integrate over the displacements. The result is 



1 



H sP = Sn,r ' S„+l, r + £ ( n ' V "> J ^ ~ r ') e( ™' r ') 



2 ^ 

n,r,r' 



+ aJ*^2 Sn.r ' S„+2,r + ^ ^(c) S„ jr • S„ ir+C , (3) 



r,c,n 



where J(r — r') is proportional to the matrix inverse of D. Now one can forget about phonons and use (Q) as a general 
model for the spin-Peierls state. 

In their low temperature phases both H s p and the antiferromagnet discussed below have order parameters. For 
the spin-Peierls model this is the staggered energy density (e(n,r)) = eo (for simplicity we assume that there is no 
antiferromagnetic order in the spin-Peierls model) . Therefore in the low temperature phase it is convenient to subtract 
from e(n, r) its average value (: e := e — eo). This leads to the following Hamiltonian 

H sP = ]T + 1 J2 : < n > r ) : J ( r - r ') : < n > r ') : + E K ^ § ^ ■ g «^+c , (4) 

r n,rr' r,c,n 

+ 2( J a + Jb)e E e(n, r) - Nel(J a + J%) , (5) 

n n n 

where J a ^ are the couplings (generated by the coupling to the phonons) to the neighbouring chains in a and b 
directions (the chain direction is chosen to be c). For defmiteness we assume in the following that J a ^ < 0. However, 
the general case can be treated by minor modifications of the formulas presented below. 

Sometimes the purely one-dimensional model (||) is used to describe the low temperature phase. We shall show 
below that such approximation always becomes very poor close to the spectral gap. 

The second model we want to discuss consists of coupled spin-1/2 antiferromagnetic chains. Their Hamiltonian is 

HaFM — J / ^(S w>r ' S„+i r) + J± S„ jr • S„ ir +a , (6) 

n,r n,r,a 

where we take J± < and a are lattice vectors in transverse directions. For simplicity we assume the transverse 
coupling to be isotropic in x and y directions. Taking the antiferromagnetic order (which we assume to be along the 
z direction) into account by a mean- field analysis of the interchain interaction one obtains the effective single-chain 
problem [[29| 

H = S„ • S n+1 -ft^(-l)^ - 2NJ ± m 2 Q , (7) 

n n 

where tuq = ((—l) n S*) and h = —4J±mQ. The remaining interchain interactions will be treated in RPA (see Section 
VI). 



II. TRANSITION TO CONTINUOUS DESCRIPTION ALONG THE CHAINS AND BOSONIZATION 

We shall discuss the spin-Peierls case of weak dimerization |(«)| « 1. In this case one can use a continuous 
description of the spin-1/2 Heisenberg chain which, for a < a c ~ 0.25 JlT]], is given by the Gaussian model. In 
the framework of this model one can express spin operators in terms of bosonic exponents (see e.g. |l3| ). Thus, for 
instance, we have 

e(x) — — — cos^^) + less singular terms . (8) 
Tra 



The value (3 = \/27r was found by Nakano and Fukuyama |14| by using the Jordan- Wigner transformation with a 
subsequent bosonization. There is a simplier way to establish this value of /?, namely, we can use the fact that the 
initial Hamiltonian has an SU(2) symmetry which must be respected by the bosonized form. Replacing u n in Eq.(|) 
by the order parameter, that is treating it as a number and using Eq.(|8|) we obtain the sine-Gordon model. Since the 
SU(2) symmetry in the sine-Gordon model is not present for general values of (3 its requirement imposes a restriction 
on the value of (3. The corresponding point in the sine-Gordon spectrum was discovered by Coleman jl2| and Haldane 
jl5| who pointed out that at f3 = \/2it there are only two breathers; the first one has the same mass as kink and 
antikink (let us call it M) and the second has the mass equal to ^/3M. Therefore at f3 — ^/2n kink, antikink and the 
first breather realize an SU(2) triplet and the second breather becomes an SU(2) singlet. 



2 



In what follows we also need the bosonized form of the spin density in the continuum (S n 
which is given by 



a S(x) 



na ), 



S(x) 

J R,,L 



1 



+ J L (x) + {- 



l) n n(a;) 
j+ - 

J R,L — 



n z (x) = sin V2tt$(x), 



2irao 
A 



exp (=Fi\/27r($ T 0) 



Tra 



i ± (x) = exp[±iV27re(a 

7ra 



(9) 



Here n a (x) are the components of the staggered magnetization, J^'^ are the currents of left and right moving fermions, 
7ia = x and A is a nonuniversal coefficient related to the bandgap for the charge excitations in the itinerant electron 
model that gives rise to the spin Hamiltonian. The field 9 is the dual of the scalar field $ and obeys d x Q(x) — H(x), 
where II is the canonical conjugate of <&. We note that (||) differ from the "usual" expressions (see p. 270-271 of [ |l3| ) 
by a shift of the bosonic field by \/ir/&. This operation interchanges sin(v / 27r < I ) ) and cos(v / 2~7r < I>), but changes neither 
derivatives of $ nor the dual field 9. 

Notice that despite the fact that e(x) has the same dimension as the z-component of the staggered magnetization, 
it is given by a different operator (sine instead of cosine). In fact e(x) is the 2fcp-component of the charge density in 
the system with a frozen charge field. As we shall see the situation is somewhat similar to that for the spin-ladder 
(see for example [fi~6|), but there are also certain subtle differences. Substituting (JsJ) into Eqs.(^) we get the following 



bosonized version of the spin-Peierls Hamiltonian H s p = J2 r -^r '' + V % 



phonon 



2 



phonon 



J dx[H 2 r + (<9 T $ r ) 2 ] + n J dx cos(V2^$ r 

a 2 r 

> / dx : cos(\Z2^$ r (x)) : J(r - r') 
n a o J 

dx K(c) ft r (x) ■ n r+c (x) , 



l - Ne 2 Q ( J a + J b ) 



(10) 



fill 



(12) 



where /i = ^^{Ja + Jb)- Note that we have kept only the most relevant terms and neglected e.g. the Umklapp term 
(for a discussion of the role of the Umklapp term see [Hi). Our strategy is now to use exact results for the sine-Gordon 
Hamiltonian (|l^) describing a single chain (see also |l8| , |l9| ) and to treat the residual interchain interaction V in the 
Random Phase Approximation (RPA). 

The Hamiltonian ( |l0| ) describing a single chain is of the form vH^ — Ne 2 ) (J a + Jb)- The spectrum of Hq is well-known 
pl| : it is described in terms of a soliton, an antisoliton, a bound state (called "breather") of mass M and a second 
breather of mass \/3M. The mass gap M of the model is due to the dimerization caused by the coupling to the 
phonons. It is related to the scale /i as follows 20 



r(j) 
*r(i) 



M 



F£(f) 
4 r(|) 



(13) 



The ground state energy density of Hq is then given by e = — tan -| , which in turn yields an expression for the 
ground state energy density of H as a function of £q- Minimization with respect to eo yields the following mean field 



expressions for eo and mass gap M. of H. 
couplings J at b and J 



(0) 



as functions of the dimensionless (nonuniversal) parameter — and the 



M=vM 



tan • 



12 

tan 



s / 2r(i) 



2r( 



A^ 
«o 



Ja + M 



| / 2r(i) \ 3 /2r(|)\ 2 A 2 |j a + j fc | 
12 \^r(l)J \ra)J a 2 



2Jk 



C 



=: C 

,\Ja + M 



\Jq + Jb\ 

2Jk 



(14) 



Here we have used v — ^JaoK for the Fermi velocity, where k is a function of the n.n.n. coupling a. The ratio of the 
constants C and C is found to be 



C_ 

C 



0.175013 . 



33V27T 

Equation makes it clear that the gap originates from the interchain interactions. 



(15) 
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III. SINE-GORDON CORRELATION FUNCTIONS AT (3 2 = 2tt 



In this section we derive exact results for various correlation functions of the sine-Gordon model (|10|) for (3 = \J2-k. 
We start by constructing a convenient basis of states for the sine-Gordon theory by means of the Zamolodchikov- 
Faddeev algebra. This is based on the knowledge of the exact spectrum and scattering matrix of the model p2|,23[. 
We then formulate the problem of calculating correlation functions in terms of formfactors and finally give explicit 
results for the first few terms in the formfactor expansion. 

The Zamolodchikov-Faddeev (ZF) algebra for the sine-Gordon model with 1 — 2ir was derived by Affleck |Tl| ] 
who suggested a representation which manifestly respects the SU(2) symmetry. As mentioned above there are three 
single particle states with mass M which form a triplet under the SU(2) symmetry. The corresponding creation and 
annihilation operators are denoted by Z+(9), Z a {9) (a = ±i,l). Here 1 denotes the breather state and ±i denote 
soliton and antisoliton states respectively. In addition there is one single-particle breather state with mass 
which transforms as a singlet under SU(2). Its creation and annihilation operators are denoted by Z^{9), Z 2 (9). As 
usual the eigenstates are parametrized by a rapidity variable 9 such that their momentum and energy are equal to 



// , = Mj sinh 9j , < , = Mj cosh , , 



(16) 



where Mj — y3M for the singlet state and Mj = M for the triplet states. By definition the ZF operators (and their 
hermitean conjugates) satisfy the following algebra 

Za(9 1 )Z b (8 2 ) = - 9 2 )Z b (9 2 )Z a (9 1 ) 

Z a {9 1 )Z 2 {9 2 ) = S a . 2 (9i - 9 2 )Z 2 {9 2 )Z a {9 1 ) 
Z 2 {9 X )Z 2 (9 2 ) = S 2 , 2 {9 X - 9 2 )B{9 2 )B(9 1 ) , (17) 

where the two-particle scattering matrices Sij (9) are given by 

sinh ( 



1 

= ± 2' 1 



S a , b (9) 
Sa.2{d) 



ism j 



sinh 9 

S Q (9 + 
( sinh 9 



- ism 

7T 



6 ' 



S (9 



S (9) 



a, b - 



1 



1 

±-, 1 

2' 



ism 7j 



\ sinh 9 — 

For the creation and annihilation operators we have 



(18) 



Z a {9 1 )Z+{9 2 ) = S {9 2 - 9 1 )Z+(9 2 )Z a (9 1 ) + 2ir5 ab 5(6 1 - 9 2 ) 

Z a {9 1 )Z+{9 2 ) = S a . 2 {9 2 - 9 1 )Z+(9 2 )Z a (9 1 ) , 
Z 2 (9 1 )Z+(9 2 ) = S 2t2 (6 2 - 9 1 )Z+(9 2 )Z 2 (9 1 ) + 2irS(9 1 - 9 2 ) . 



a,b=±-,l 
2 



(19) 



From ( |lq ) it follows that S'i ! i(0) = — 1 and Si j(oo) = +1. Therefore particles with close momenta behave like free 
fermions and particles far apart in momentum space behave like free bosons. 

We note that the soliton S-matrix So{9) has simple poles at 9 = and 9 = i^-. These poles correspond to the 
two breather bound states. In other words a soliton and an antisoliton can form a bound state of either mass M or 
mass y/SM. Although the S-matrix of light breathers is the same only the pole 9 = z-| corresponds to a bound state 
- the heavy breather. The pole at 9 = i^- is redundant [Q. The soliton- breather S-matrices S a .i(9) and S a>2 (9) 
exhibit soliton poles at 8 = ± and 9 = ± i-| respectively. All other poles are redundant and do not indicate 
the presence of bound states. ^From the analytic properties of the S-matrices ( |T8| ) we deduce relations between the 
ZF operators. For example we find 



Z 2 {- 



lim 



£i(0 2 )2i(0i) 



(20) 



Such relations between the ZF operators play an important role in what follows. 

States in the Fock space are constructed by acting with the operators Z\(G) on the vacuum state |0) 
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\e n ...6 1 ) £n ... ei =Zl(6 n )...Zl i (9 1 )\0), (21) 

where Sj = ±g, 1, 2. We note that ( |l7| ) together with ( pl| ) implies that states with different ordering of two rapidities 
and indices £i are related by multiplication with 2-particle S-matrices 

\6 n . .. 8 k 9 k+ i . . . 6x)e n ...£ k e h+l ... ei = Se k ,e k+1 (Ok ~ k +l)\0 n . . . 9 k+1 6 k . . . 9l)e n ... ek+1 e k ...e 1 (22) 

The resolution of the identity is given by 

1 = EE/ ■ (23) 

n=o Ej J y ) ■ 

The formfactor approach is based on the idea of inserting (|23|) between the operators in a correlation function 



(O(*,*)Ot(0,0)> = f:5: / ^^exp (i± PjX - 

n=0 e; V ' \ 7=1 



jt K0|O(0,0)|^...^) en ... ei | 2 , (24) 



and then determining the formfactors 

F°(9 1 ...d n ) £l ... En :=(0\O(0,0)\6 n ...e 1 ) en ... ei (25) 

by taking advantage of their known analytic properties. 

From a physical point of view we are interested in the Fourier transforms of the connected retarded 2-point corre- 
lators of cos ^/2tt^ and sin-v/^^. Their formfactor expansions are of the form 



/oo />oo 
dx / dt e' i( " +le)t ~ l ^([cosV2^$(M),cos\/2^$(0,0)]) 
-oo JO 



(27r)"n! 1 v ' r ^ ei - s »i w-]T,M, cosh0, -He w + M, cosh 0,- + ie ' 

(26) 

Here we have reinserted the Fermi velocity v and lattice spacing ao- The Fourier transform D sm (ui, q) of the connected 
retarded 2-point correlator of sin-\/27r3> is the dynamical staggered susceptibility and will also be denoted by x"(w, q). 

In order to implement the formfactor expansion it is very useful to note that (like for general values of 0) operators 
from different representations behave differently under the charge conjugation transformation 

C$C _1 = -$ , 
CZ ± i (e)C- 1 = Z±i(6), CZi{6)C~ l = -Zi{0) , 

CZ 2 {0)C- x = Z 2 (6) . (27) 
These transformation properties imply the following expansion 

sin[v&(T,a:)]|0) = F 1 f ^e- iM ( tcoshe - xsinh ^(0)|O) 

J 27T 

j^l ^2 e ~iM[t (cosh 0i+cosh 2 )-x(sinh 9i+sinh 9a)] x 

2tt 2tt 

x U(9 ll d 2 )[Z+(0 1 )Z+ 1 (0 2 ) - Z+ 1 (d 1 )Zt(6 2 )} + ... (28) 

2 2 2 2 

where U(8i,9 2 ) — U(9 2 ,8i). Significantly, due to the SU(2) symmetry transverse components of the staggered 
magnetization have the same correlation functions as n z stag (this is clear from the SU(2) symmetry of the Hamiltonian 

(|J)). Consequently we conclude that Sine-Gordon 2-point correlation functions (for (3 = \/2tt) of sin\/27r# and 
cos-v/^^ are the same as 2-point correlators of sin\/27r$. 

The current operator d x <f> is also odd in $ and therefore its expansion must begin with Zq such that at small q we 
have 
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((g( W) g)g(-w,-g)>)~— _p— -— + ... (29) 
uu^ — — M z 

a o 

where dots denote terms which have nonzero imaginary parts at higher energies. 

As we will see the threshold of the dynamical spin susceptibility is equal to M for both q — and q = n. This is a 
distinct feature of the alternating chain. It is related to the fact that kink and antikink create a bound state of the 
same mass. Recall that for the ladder chain (or S = 1 antiferromagnet for this matter) where particles do not have 
bound states, the value of the energy threshold at q = is twice that at q — n (see Ref. fHI ). 

At frequencies smaller than (1 + \/S)M the only contributions to the imaginary part of the magnetic susceptibility 
come from the first breather and kink-antikink pairs. Kink-antikink form factors can be calculated in the Sine-Gordon 
model (for any value of the coupling (3) along the following lines |^,^6|. Let us denote by S + {8) (S_(8)) the S- matrix 
eigenvalue corresponding to positive (negative) C-parity obtained by diagonalising the kink-antikink scattering: 

sinh &(9 + in) cosh £ (8 + in) 

;S (8), S-= . 5 (g) , (30) 



sinh || (9 — in) ' cosh ^ {9 — in) 



where £ = 87T Z.^ and 



f dx sinfasmli^i 

S {9) = -exp(-i / , \ ) . (31) 

Jo x C osh ^ sum ^ 

Then, general unitarity and crossing arguments imply that the corresponding kink-antikink form factors F±{8) are 
solutions of the following system of functional equations 

F±(-9) = S ± (9)F ± (9) (32) 
F± {9 - 2in) = ±F±(-8) . (33) 

The "minimal" solutions of these equations are 

F + (9) = . *fl 9 . , F (9) 
smh(9 + m) 

F^O) = -J^-- F (9) , (34) 



cosh(0 + i7r) 



■2S. 



, da;smhf(l-|)sm 2 i^^ x , ^ 

F Q (9) = sinh - exp( / . 2 \ f *' t _ ^ ) . (35) 



where Fq is given by 



2 Jo x sinh §£ cosh f sinh x 

By minimal solution we mean a solution containing only the expected bound state poles in the physical strip and 
with the mildest asymptotic behaviour at infinity. This prescription determines the minimal solution uniquely. An 
infinite number of non-minimal solutions corresponding to all operators in the theory which are local with respect to 
the solitons are obtained multiplying the minimal solution by an analytic function of cosh#. However, if we require 
the form factor to be power bounded in the momenta and to have only the bound state poles, we conclude that we can 
actually multiply the minimal solution only by a polynomial in cosh#. For a given operator, it is possible to put strong 
constraints on the asymptotic behaviour of its form factors, and then on the degree of the allowed polynomial p7| . 
In the sine-Gordon model this procedure is complicated by a non-trivial behaviour of correlators in the ultraviolet 
limit. Nevertheless, the result is that for the operators cos and sin the allowed polynomial is of the zero degree, which 
means that their form factors coincide with the minimal ones. The same conlusion can be reached in a simpler way 
going to the free fermion point £ = n, where the form factors of sin and cos can be easily computed remembering that 

cos/3$ - £^3"$ - Jf, , (36) 

and that the sin is related to the elementary field by the equation of motion. 

For the operators cos (3$ and sin/3$, at the specific value of the coupling we are interested in, we find (812 = 81 — 82) 



G 



|sinV2^iA)- + = ^(^^ ZhwJ/^ 6 ^ = FSm{6) - + (37) 
2^\0 1 ,e 2 )-+ = iV3(2d)Z^ '■ l^ll 2 ^ 2 C{0i2) = F c ° s (9 12 )_ + (38) 



C(6») = csinh0/2exp da; 



sin 2 [x(9 + in)/ 2] cosh7r:r/6 
x sinh nx cosh 7rx /2 



, 1 Z" 50 sinh x/2 cosh x/6 ] , 3 

c= (12)3 exp<^ / d.x 2 , \~ 3.494607 , d = — w 0.136629 , (39) 

2 Jo xcosh a;/2 J 2nc 



where the relative normalisation between the two operators can be fixed exploiting the asymptotic factorisation of 
form factors discussed in Ref. ]2S|l . We note that £(#) is to be analytically continued using the relation 

C(e)s (e) = c(-e) . (40) 

The additional factors d and c in (|J7|) and ([58]) have been introduced in order to simplify the reduction of multiparticle 
formfactors using the annihilation-pole condition (for soliton formfactors) 



Ota a_ ^ . . £- £ 2™-i 

£2n 



— F (#1 ■ • ■ &2n-2) e ' 1 ...e' 2n _ 2 S l 

62n—62n-l=iK 



iResF°(e 1 ...9 2n ) £l .. 

X jljgj . . . 5 E l™Zl — St^e^ 1 ' 1 (02n-l — ) ■ ■ • S e lZ-l',elZ-l (02n-l — 02n-2)j i (41) 

where S|'f (0) = S"/ (0) = Si and all other components are zero. Multiparticle formfactors are discussed in some 
detail in Appendix A. 

The form factor ( |37j ) has a pole at #12 = — 2i7r/3 corresponding to formation of a bound state - the first breather. 
The breather formfactor F\ is given by the residue of (STi) divided by the three-particle coupling: 



,0 32 / f°° dx sinh 7rizi/6 sinh 7rx/3\ 

LFi 2 = — exp -2 / — '- — tz-)Zks 0.0533Z . 42 

87H \ Jo xsmh7ra;cosh7rx/2 / 

Similarly ( p8[ ) has a pole at 6> 12 = — i7r/3 corresponding to the second breather. The absolute square \F 2 | 2 of the 
breather formfactor is found to be 

l9 3^ / f 00 dx cosh ttx/6 sinh 2 7rir/3\ . , 

F 2 2 = — r exp -4 / — ' j-J— Z w 0.0262Z . (43) 

87r^ \ Jo a; sinh to cosh irx/ 2 / 



IV. DYNAMICAL SUSCEPTIBILITIES FOR A SINGLE ALTERNATING CHAIN 

The expression for the imaginary part of the dynamical staggered susceptibility ~x"(uo,q) at s 2 = uj 2 — ^q 2 < 
(1 + V3) 2 M 2 is given by 

Zm X "(w,q) = 2 7 r|F 1 | 2 ,5( S 2 - M 2 ) + 2SRe |Jr 7f (S !!t: 2 |2 . ( 44 ) 



where 0(s) = 21n(s/2M + s 2 /AM 2 — 1). Note that all other formfactors do not contribute to this expression in 
the specified range of s as their thresholds are above (1 + \/~3)M . Also the normalization Z enters ( ^4]) only as an 
overall factor. Since the function ((8) vanishes at 9 = 0, the entire formfactor is also finite. Thus the two-particle 
contribution to <?) exhibits a square-root singularity at the threshold as a function of s. 

The breather and ss contributions to the real part are found to be 

2|^i| 2 „ f°° d0 s 2 - AM 2 cosh 2 f 



3?e X "(w,g) = -Me- — ^JJ 2 / ^-—^ |F s,n (0) + _| 2 , (45) 

AV ' S 2 -M 2 + IS J 7T ( S 2_4 M 2 cosh 2| )2 + e 2l V ' 

where the factor of 2 stems from the sum over + and — . In Fig. |l| we plot both the imaginary and real parts of x" 
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FIG. 1. Imaginary and real parts (in units of jj?) of the dynamical staggered susceptibility as functions of s = 
for q ~ n. The dashed line depicts the single-mode approximation that takes into account only the first breather. 

It is straightforward to repeat the above analysis for the current operator d x <!> using the explicit expressions for 
the formfactors given in p6| . The contribution of the first breather leads to to (^9|) with some normalization factor. 
As there is very little spectral weight at q ss we concentrate on q w ir and do not repeat the above analysis for the 
current operator. 

For practical purposes it is convenient to have an expression interpolating between the small q ( p9|) and q sa tt ( |44| ) 
behaviour. Such an expression giving the dynamical spin susceptibility in the entire range of q at frequencies below 
the continuum may look like 



g{q) sin 2 (<7/2) 



(46) 



w 2 -4sin 2 g-M 2 "' ' 

a o 

where g(q) is a smooth function interpolating between the normalizations at q — and q = tt. The mode lu = 

A /t- sin 2 q + M 2 is separated from the particle continuum by the gap of order of M, 
V a o 

Let us now turn to the two-point correlator of cosines. The contributions of the second breather and the soliton- 
antisoliton continuum are given by 



%mD cos (Lu,q) = 2tt\F 2 \ 2 8{s 2 - 3M 2 ) + 25ic 
where the ratio of the single particle residues is universal: 

sinhx/3 



F cos [8(s)l 



sys 2 



AM 2 



7 : 



exp 



dx- 



x cosh 2 | 



0.49131 



(47) 



(48) 



Note that the threshold of the breather-breather continuum is also at s = 2M. The corresponding contribution is 
taken into account in Appendix A. The analogous contributions to the real part of D cos are given by 



MeD cos (Lu,q) = -5Rc- 



3AF 



is 



de 



AM 2 cosh 2 



tt ( s 2_4M 2 cosh 2 f) 2 +£ 2 



|f««(0) + _|2 



(49) 



The remaining integrals in ([14|),(|45|),(|47|) and (|49|) have to be calculated numerically. We find that at small s the 
contributions of the two-particle continua to the real parts of both correlators are of the same magnitude as the single- 
particle contributions from the breather states. As far as a single chain is concerned a single-mode approximation 
taking into account only the one-particle states is therefore very poor at small s. 
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FIG. 2. Imaginary and real parts of D cos (in units of -^) as functions of s = , u 2 — ^q 2 for q ~ tt. 



V. RPA ANALYSIS OF THE INTERCHAIN INTERACTIONS 

Let us now take into account the interchain interactions (both of spin and staggered energy densities) in (|]). This 
is accomplished through an RPA analysis along the lines of |30| ■ RPA becomes exact in the limit of an infinite number 
of neighbouring chains. 

In the RPA we obtain the following expression for the correlation function of energy densities 



Xe(s, k) = ((e(-uj, -q; -k)e(ui, q; k)}) = 



D cos (s) 



1 - D cos {s) J(k) 



where J(k) — 2[\J a \ cos(k x ) + | J&| cos(k y )]. Similarly the dynamical staggered susceptibility is given by 



X zz (s, k) = «S z (-c, -q; -k)S z {u, q; k))) = 



D sm {s) 



1 + D sm (s) K{k) 



(50) 



(51) 



where K(k) = 2[K a cos(k x ) + K\> cos(fc y )]. Note that in the present approximation x e (s,fc) is only affected by the 

interchain interactions of staggered energy densities whereas x zz { s ik) only "sees" the interchain interactions of spin 
densities. The reason for this decoupling is that in the sine-Gordon theory describing the individual chains 



(e(t,x)S* (0,0)) = 



(52) 



because e(t, x) is even under charge conjugation whereas S" z (0, 0) is odd. By rotational invariance this implies 
(e(t, x)S(0, 0)) = 0. Note that the RPA is particularly simple as we have taken into account only interchain in- 
teractions of the staggered part of the spin density and neglected the smooth part as being less relevant. If we take 
these subleading terms into account the RPA acquires a matrix structure like in p9| as the sectors q ~ tt and q ~ 
become coupled. An RPA analysis then requires the calculation of formfactors of the current operator. We will discuss 
this refined RPA in a separate publication M]. 

The response functions x aa ( s i k) and Xe( s J k) can be easily calculated numerically for given values of J a ,b and K a ^ 
by using the expressions for D sm (s) and D cos (s) obtained in section IV. We note that the pole in the dynamical 
magnetic susceptibilities x aa ( s i k) corresponds to the light breather which has quantum number S z = 1, whereas the 
pole in Xe( s i k) is due to the heavy breather with S z = 0. The magnetic mode can be measured directly by neutron 
scattering whereas the S z — excitation can be probed by measuring the phonon spectrum which will exhibit a 
softening. 

In order to visualize our results we now plot them for a particular choice of parameters. Being aware that our 
theory probably cannot be applied to CuGeC>3 where a > a c we nevertheless choose these parameters to reproduce 
the dispersions of magnetic excitations in that material. We take n ~ 0.8, M « A.58meV and ^jpr- ~ 3.08, 







1 Tp w 0.25. The value of k is chosen such that the dispersion in z-direction reproduces the experimental fit of 
j|, whereas the other conditions follow from the experimental band-gaps for fc-vectors (0, 1, |) (~ 2meV), (0,0, |) 
(« 5.7meV) and (§, 1, \) (« 2.6meV). 

Because of Lorentz-invariance the energy u> and the z-component of the momentum q only enter in the combination 

s = y^j 2 — ^q 2 - In Fig. [| (a) we plot the spin- wave dispersion in ^-direction (k y = 0, k x S (0, 7r)) and in Fig. [3] (b) 
in rr-direction (k x — 0, k y £ (0,7r)). 
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FIG. 3. Spin-wave dispersions in x (a) and y (b) directions. 



We see that the single-mode approximation (SMA) in which all multiparticle contributions to the dynamical sus- 
ceptibilities are neglected gives essentially the same result as the exact treatment. We note that (by construction) the 
fits Q to the experimental results are essentially identical to the SMA as far as dispersion relations are concerned. Let 
us now turn to the multiparticle continuum. The imaginary part of the dynamical staggered susceptibility is directly 
measurable by Neutron scattering. The position of its poles yields the dispersion discussed above. The incoherent 
part (as a function of s and k XjV ) is plotted in Fig. || (a) and (b) respectively. 



Imx z 




Im % 




0.05 



(a) (b) 

FIG. 4. Imaginary part (in arbitrary units) of \ zz (u,q; k) in x (a) and y (b) directions. 
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In the RPA the 2-particle continuum starts at s = 2M w 9meV. This is in disagreement with experiment for 
CuGe03. This may be because a > a c or because the dispersion in ^-direction is rather strong, which in turn forces 
M to be large in order to get a reasonable fit to the experimental data. 



VI. ANTIFERROMAGNETIC ORDER 



In the continuum limit the model (pi) is also equivalent to the sine-Gordon model with (3 2 = 2n 29 1. The bosonization 
formulas are now 

S(x)=3 R {x)+3 L (x) + {-l) n H(x) , 

7= (d x ® t n) , J+ L = J— CX p ( T tV2^^ T e) 
'2n ' Zirao V 



J R,L — 



n z (x) = — cos V2tt<&(x), ii ± (x) = — — exp[±i\/27r6(x) 

7rao 7rao 



(53) 



Since the staggered magnetization is proportional to sin(\/27r$) the bosonized single-chain Hamiltonian is given by 



H 



dx [TT 2 (x) + (d x <S>{x)f 



hX 
Tra 



dx sin V27r$ - 2NJ_ L m„ 



(54) 



Following through the same steps as in the spin-Peierls case we find the following mean-field results for staggered 
magnetization and mass gap 



m 



J± 



J 



M = C'\J ± \ , 



(55) 



where the ratio of C and C is given by (|15|). This relation can be used to determine the transverse coupling J± in 
terms of J and the directly measurable quantities M and too as follows. The gap M is equal to u>(tt,0,tt) which 
is found experimentally to be 11.0 ± O.bmeV in KCUF3 |^|. The average magnetic moment is too ~ 0.27 B and 
J ~ 53.17 ± 0.25meV 0. Using these values we find 



\J±\ 



1 / CM 



J V C'toq 



0.96meV 



(56) 



Let us now turn to the correlation functions for a single chain. The correlator ((S Z S Z )) at small q is still given by 
Eq. (p9|) , but for q w it is given by Eq.(|47|). Therefore around q = w the pole is at s = \/3M 



X **{w,q) = D c ° s {s) 



2Fi 



3M 2 + £(TT- q y 



incoherent 



(57) 



where "incoherent" denotes multiparticle contributions. Some of these are determined exactly in Appendix A. A plot 
of D cos (s) is shown in Figj|. 

Correlation functions of transverse components of the staggered magnetization are given by Eq.(^4|) and have a 
pole at s — M 



X xx (u,q)^D sm (s) 



2Ff 



M 2 + £(n-q)2-LjZ 



incoherent 



(58) 



Here we have used the fact that in the continuum model the correlation functions of cos y'27rd, sin \/2tt6 and sin v 27r$ 
are equal due to the SU(2) symmetry present at = y/2n. A plot of D sm (s) is shown in Fig. 0. 

The difference in the correlation functions (|57|) and ( ^8[ ) is obviously related to the broken rotational symmetry 
of the Hamiltonian (0). Next we take into account the interchain interaction by an RPA analysis. This yields the 
following expression for the longitudinal dynamical susceptibility 



X zz (uJ,q,k) = 
X xx (uJ,q,k) = 



D cos (s) 



1 - 2|J_l|(cos k x + cos k y )D cos (s) 
D sm (s) 

1 - 2| Jj_ I (cos k x + cos k y )D sm (s) 



(59) 
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Here we again have taken into account only the staggered part of the spin density as it gives the most relevant 
contribution to the interchain interaction. As a result the RPA expression for the susceptibilities are of scalar rather 
than matrix form [p9[ . The transverse susceptibility must have a pole at the Neel wave vector (0,0, n) as the spin 
SU(2) symmetry is spontaneously broken. This leads to the requirement that 

DSm(0) = ^"°- 125 ° 9 ^' (60) 

which fixes the normalization Z in terms of the transverse coupling and the breather mass as Z « 1.999^. The 
normalization of the correlator of cosines then follows to be 

^•(0)-^. (61) 

The Goldstone mode associated with the zero energy pole in \^ is a spin-wave moving in z-direction and its 
dispersion is found from the singularities of \ . Due to Lorentz invariance D sm depends only on s rather than on u> 
and q independently. This immediately implies that the spin-wave dispersion for q w tt is 

Jfa %) = 4(* - if + M 2 (l - C ° Sfc - +C0S S « ^(vr - qf + M\l CQSfc - +C0S ^) . (62) 
Gq 2 4 2 

This is in very good agreement with experiment being almost identical to the fit used in |^] . In Fig. [5] (a) we plot the 
spin-wave dispersion for k y = and k x £ (0, n). 




0.0 1.0 2.0 3.0 4.0 2 



(a) x (b) 

FIG. 5. (a) Spin-wave dispersion as a function of k x for k y = 0. (b) Imaginary part (in arbitrary units) of x ± for s > 2 as a 
function of k x for k y = 0. 

We see that the SMA works extremely well for all values of k x . The imaginary part of the dynamical susceptibility 
is directly measurable by neutron scattering. We find 

c* ±/ f\ n x( s2 1 , cos fc^ + cos fc. y . 

3mx (w, g, k) = -—- 8{j^ - 1 H + incoherent , (63) 

where we have used the SMA to get the delta- function part. The incoherent part is plotted in Fig. || (b). We see that 
there is in general no singularity at the threshold of the light breather-heavy breather continuum except at k x — > tt 
where Sm^ 1 ^, q, k) = 3mD s ™(u, q) so that we recover the pure 1-D result. The situation for the soliton-antisoliton 
continuum is analogous. 
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Let us now turn to the longitudial susceptibility. In the SMA there is a pole in \ zz a ^ 

2 

L0 2 {q, k) = ^(n - q) 2 + M 2 [3- ^- (cos k x + cos k y )} , (64) 
a 2 

where 7 is given by (pL8|). This is compared to the exact result for the case where k y = in Fig. || (a). We see that 
the corrections to the SMA result are very small. 



1.70 



1.65 



1.60 



1.55 



Imx 



Exact 
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0.0 



1.0 



2.0 
k 



3.0 



4.0 



s/M 



FIG. 6. (a) Dispersion of the longitudinal mode as a function of k x for k y = 0. (b) Imaginary part (in arbitrary units) of 
X zz for s > 2 as a function of k x for k y — 0. 

Using the SMA (which we know from Fig. ^(a) to be an excellent approximation) to extract the coherent delta- 
function part we find 



-* 7T"7 S 7 

3mx zz (w,g, fc) = — - r d(jj^2 ~ 3 + 2 ( C0S Mr + cosk y )) + incoherent . 
The incoherent part is plotted in Fig. §(b). 
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APPENDIX A: MULTIPARTICLE FORMFACTORS 



(65) 



In this appendix we consider multiparticlc formfactors. We start with 2-soliton 2-antisoliton formfactors and 
explicitly derive the related three and two-particle formfactors. The extension to n-soliton n-antisoliton formfactors 
(n = 3,4,...) is straightforward and will not be discussed here. The formfactor expansion for 2-point correlation 
functions is found to be rapidly converging and for small s it is essentially sufficient to take into account 2-particle 
formfactors only. Note that most of the formulas below are to be understood in terms of analytic continuation of 
A useful formula is 



\((6 + ia)C(0 - ia)\ 2 = j(cosh( 



cos a) exp < 4 



coshx/6 , x6 , f-a . 

ax — — r^r(l — cos — cosh x) 

x smh x cosh | ir ir 



sinh 2 8 cos 2 a + (cosh 8 sin a — sin -|) 2 
sinh 2 8 cos 2 a + (cosh 8 sin a + sin -|) 2 



(66) 
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Multiparticle formfactors of (quasi)local operators were studied in detail by Smirnov pq] . From his work the 
2-soliton 2-antisoliton formfactor for cos \/2ir& is straightforwardly extracted 

F cos {6x, 62,03, 64) — ++ = (O|cosv&|0i,03,03A)— ++ = 

2.(2d)Vz 3 (X>^ fE^)E e9i+fl5 'II^-^K IE/i >< 

\i = l / \m=l / i<j i<fc 

x cosh^3 + ^-^^ 2 )firi sinh3 ^ (67) 

i— 1 j= 3 J 

Orderings other than h+ are obtained from (^) by using the generalization of ( |33| ) e.g. 

F cos (6i, 82, 6*3, $4) — 1 — h = £o(#3 — #2) F cos (8i, 63, 62,04) h+ . (68) 



It is easy to verify that the soliton-antisoliton formfactor of cosv27r$ is obtained from (67) via the annihilation pole 
condition ([ll]) . The 2-soliton 2-antisoliton formfactor for sin \/2tt$> is very similar to the one for cos \/2tt$> 

F sm (6 1 ,8 2 ,6 3 , 4 )— ++ = -itanh |(d 3 +6 i -0 l - 6 2 )F cos (d 1 ,d2, 63,64)—++ ■ ( 69 ) 



The residue at the annihilation pole (times i) now yields the soliton-antisoliton formfactor of sin V27r$. 

Breather formfactors are obtained from the residues of (|6^) at its poles. In the soliton-antisoliton-even breather 
sector we find 



F cos (8 1 ,8 2 , 63) < 2 = -27ri-f4 =f C(6 12 )C(-i-)C (6 13 + i- 

1 ' ' 37 + 2i3l cosh ^ cosh 3031 cosh 30 32 3 AV 13 6 ; 

x C(*13 - iJ)C(^3 + iJ)C(^23 - i£)[e fll + e 02 + V3e^] 
6 6 6 

x [ e -^+e-^+V3e- e3 ][e 6 ^ + V3e^(e^+e^) + e 2e -]. (70) 
The corresponding formfactor for sin \/2ir& is 



F sm {6 l ,6 2l 63)-+2 = -icoth^i F c "(0 1 ,0 a ,0 3 )_ +a , (71) 
and different orderings are obtained by the appropriate generalization of d33| ) e.g. 

F cos {6\, 62, 63)2 — h = F cos (02, O3, 81) — V 2 Si, 2(^21) 5*1,2(031) ■ (72) 

The residue at the annihilation pole (times i) in F cos (8i, 62, 63)2 y gives the heavy breather formfactor F2. The 

corresponding sin formfactor has no annihilation poles. In the soliton-antisoliton-odd breather sector we obtain 

(2d) 2 J~Z e -i?( i+«2)-03 2tt tt 

F^8 M ^ +1 = - 2 , ^| sinh 4 sinh 3, 3ismh3g23 CC^)C(-if)C(^s 

x C(013 - i|)C(^23 + i|)C(^23 - + e e2 + e 6 *} 



x [e _ei + e~ e2 + e~ e3 ][e ei+e2 + e Bi {e ei + e" 2 ) + er 



20 ; 



= icoth^i F* n (9 1 ,6 a ,0 3 )-+i . (73) 



From the residues at the poles of ( |73| ) we can derive the breather-breather formfactors: 

A{2d) 2 VZ (cosh ^f) 2 [cosh <9 12 + |] 



F cos f0i,# 2 )ii = -2tt- 



3 1 (sinh30i 2 ) 



x C 2 ^2i)C 2 (-iy)C(^i + iy)C^2i - iy)5 o (02i + i^)S (6 21 - ig ) • ( 74 ) 
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This is identical to the soliton-antisoliton formfactor F 00 * (61,62) \- as can be proved by direct calculation. Some 

sinh 36 
sinh sin 



useful identities are S (6 + if )S (6 - if) = S {0), 2C(0)C(# - »*■) = jff' and 



, , dx cosh f - cosh § \ 2^/2 , . 

-— 7 s -— r^H— ■ (75) 



r sinh x cosh ^ 



The other breather-breather formfactors are given by 



i(2d) 2 VZ (cosh ^f) 2 [cosh 12 + |] 



F si "(^,^)2i = -27r 



(cosh30i 2 ) 2 

xC 2 (MC 2 (-i|)C(^i2 + i|)C(^i2-i|) , 

4(2d) 2 VZ (cosh0 12 + a^)[l + ^ CO sh(9 12 ] 



(cosh30i 2 ) 5 



x C(-i^)C(-iy)C(«i2 + i g )C(»i2 - i|)C(^i2 + i^)C(^i2 - ig ) . (76) 

The special values of C at the breather poles are given by C(~if ) ~ -1.10184i and CMif) ~ -2.72272L We note 
that C(-if )C(-i^ L ) = -3- 

It is apparent that n-particle formfactors depend only onn-1 independent rapidity variables. This fact can be 
used to essentially simplify expressions like (pq) for correlation functions. For example, 2-particle contributions for 
lu > are given by 



D cos (co,q) 







2 — part. 








2 — part. 





\F cos (6) 2 



s 2 _ AM 2 cosh ^ | + fe s 2 _ 12 m 2 cosli I + is 



2|F sin (6»)+_| 2 2|F sm (6>) 



21 



4M 2 cosh 2 f + s 2 _ 4M 2 (1 + & cosh 0) + ie 



(77) 



where we also have made use of various symmetry properties of the formfactors in order to perform the sum over ej . 
Similarly the contribution of formfactors involving one soliton, one antisoliton and one (light) breather of type 1 can 
be brought to the form 



00 66 f°° d6 



OO 



- cos (Mi2)-+i| 2 , (78) 



2ttJ_ 00 2tt s 2 - M 2 (1 + 4 cosh cosh 2f + [2cosh2f] 2 ) 



where 



P cos mf) , 9 (2d) 2 (2 cosh e -f + e e )(2 cosh e -f + e -«)(2 cosh 6 -f + 2 cosh0) 2tt 
F (0, 12 )_ +1 - -2*-^ S inh3^sinh3(0+^)sinh3(0-^) « 

x((|-«- if )C(^ - 6 + if ) C (^ + - if )C(^ + 6 + if). (79) 
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